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ABSTRACT 


A  recently  developed  method  for  predicting  propeller-hull  interaction,  in 
which  a  numerical  method  for  calculating  the  viscous  flow  over  the  stern  and 
in  the  wake  of  a  ship  is  coupled  with  a  propeller-performance  program  in  an 
interactive  and  iterative  manner  to  predict  the  combined  flow  field,  is  vali¬ 
dated  by  performing  comparisons  between  computational  results  and  extensive 
available  experimental  data  for  propeller-shaft  configurations.  The  steady- 
flow  results  are  in  excellent  agreement  with  the  data  and  show  that  the  pre¬ 
sent  procedures  are  able  to  accurately  predict  many  details  of  the  flow 
field.  The  dependence  of  the  flow  field  on  propeller  loading,  including  the 
formation  of  the  hub  vortex,  and  the  influence  of  hub  length,  are  accurately 
simulated.  Also,  the  robustness  of  the  solution  procedure  is  demonstrated  by 
performing  calculations  for  off-design  (large-loading)  conditions.  The 
unsteady-flow  calculations,  which  simulate  the  fanning  action  of  a  rotating 
finite-bladed  propeller  and  are  in  reasonable  agreement  with  the  experimental 
data,  point  out  the  difficulties  of  accurately  resolving  the  complex  blade-to- 
blade  flow  and  the  need  for  investigating  alternative  approaches. 
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I.  INTRODUCTION 


For  practical  reasons  marine  propellers  are  located  at  the  stern  of  a 
ship;  consequently,  they  operate  in  the  thick  stern  boundary-layer  and  near 
wake.  The  flow  field  around  the  propeller— hull  combination  is  unsteady,  three 
dimensional,  and  turbulent.  It  is  also  interactive,  insofar  as  the  propeller- 
induced  flow  is  dependent  on  the  hull  flow  which  is  itself  altered  by  the 
presence  of  the  propeller.  Historically,  methods  for  the  analysis  of  the 
propeller-induced  flow  and  those  for  the  hull  boundary  layer  and  wake  have 
been  developed  separately  with  relatively  little  attention  given  to  the  inter¬ 
action. 

The  propeller-induced  flow  is  usually  calculated  under  the  assumption 
that  the  propeller  operates  in  an  infinite  ideal  fluid,  but  with  a  specified 
spatially  varying  potential  inflow  which  represents  the  hull  boundary  layer 
and  wake.  Quite  advanced  inviscid-f low  methods  have  been  developed  to  imple¬ 
ment  this  propeller  theory.  Consistent  with  the  underlying  assumptions  of  the 
theory,  these  methods  have  been  found  to  perform  best  for  uniform  inflow 
conditions.  A  complete  evaluation  of  the  theory  is  made  difficult  by  the  lack 
of  knowledge  of  the  effective  inflow  which  is  usually  assumed  to  be  the  nom¬ 
inal  wake  of  the  bare  hull.  In  the  absence  of  reliable  methods  for  the  pre¬ 
diction  of  the  latter,  recourse  is  had  to  model  experiments  in  towing  tanks. 
The  flow  over  the  stern  and  in  the  wake  cannot  be  treated  by  the  usual  thin¬ 
boundary-layer  equations  and  is  influenced  by  both  appendages  and  free-surface 
waves.  The  theoretical  developments  on  this  topic  have  generally  focused  on 
the  bare  hull  and  without  the  effects  of  the  free-surface  waves,  and  it  is 
only  very  recently  that  the  level  of  sophistication  required  for  handling 
practical  bare  hull  geometries  has  been  reached. 

Recently,  the  authors  (Stern  et  al.,  1985  and  1986),  have  taken  a  com¬ 
prehensive  approach  in  which  methods  of  modern  computational  fluid  dynamics 
have  been  brought  to  bear  on  the  problem  of  propeller-hull  interaction.  The 
partially— parabolic  method  of  Chen  and  Patel  (1985)  for  calculating  ship— stern 
flow  has  been  coupled  with  a  propeller-performance  program  in  an  interactive 
and  iterative  manner  to  predict  the  combined  flow  field.  A  body-force  distri¬ 
bution  is  used  to  represent  the  propeller  in  the  partially-parabolic  method. 
The  overall  computational  method  is  applicable  to  the  most  general  situation 
of  unsteady  turbulent  flow  about  propeller-driven  three-dimensional  bodies. 
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Some  calculations  for  axisymmetric  and  three-dimensional  bodies  were  reported 
in  Stern  et  al.  (1985,  1986).  The  latter  reference  also  provided  typical 
results  for  the  simple  geometry  of  a  propeller-shaft  configuration. 

This  report  is  concerned  specifically  with  computations  for  propeller- 
shaft  configurations.  In  this  case,  the  upstream  body  is  simply  the  propeller 
shaft.  Although  this  is  clearly  an  idealization  of  the  practical  circum¬ 
stance,  it  will  be  shown  that  the  flow  field  exhibits  all  the  distinctive 
features  of  interest.  Furthermore,  the  simplicity  of  the  geometry  makes  it 
possible  to  identify  and  isolate  the  important  features  of  the  propeller- 
induced  flow  field.  No  doubt,  it  is  for  this  reason  that  the  most  extensive 
experimental  data  that  is  available  for  propeller  flow  fields  is  for  such  a 
configuration. 


II.  COMPUTATIONAL  METHOD 

Consider  the  boundary  layer  and  wake  of  a  ship  which  is  advancing  at 
constant  speed  UQ  under  the  action  of  a  propeller  rotating  at  constant  angular 
velocity  to  .  The  free-surface  waves  are  ignored  and  the  free  surface  is 
considered  a  plane  of  symmetry  (zero-Froude-nuraber  approximation).  The  flow 
field  may  be  divided  into  four  regions,  as  depicted  in  figure  1.  In  region  I, 
which  is  outside  the  boundary  layer  and  wake,  the  flow  is  considered  inviscid 
and  irrotational .  Region  II  contains  the  propeller.  Region  III  is  the  thin 
boundary  layer  on  the  hull,  and  ends  at  a  station  x*  where  the  boundary— layer 
approximations  are  no  longer  valid.  Region  IV  (x  >  x*)  includes  the  thick 
boundary— layer  and  wake.  Herein,  the  primary  focus  is  on  region  IV  and  the 
interactive  analysis  to  obtain  the  combined  propeller-hull  flow  field. 

A.  Method  for  Calculating  Ship— Stern  Flow. 

In  the  thick  boundary  layer  over  the  stern  and  in  the  wake  the  usual  thin 
boundary-layer  assumptions  are  no  longer  valid  and  the  review  of  Patel  (1982) 
indicates  that  the  partially-parabolic  Reynolds  equations,  in  which  only  the 
streamwise  diffusion  terms  are  neglected,  are  more  appropriate.  The  primary 
effect  of  an  operating  propeller  on  this  flow  field  is  to  accelerate  the  fluid 
in  the  vicinity  of  the  propeller  in  both  the  axial  and  tangential  directions, 
and  this  reinforces  the  partially-parabolic  assumptions.  Recently,  Chen  and 
Patel  (1985)  have  developed  a  method  for  solving  the  partially-parabolic 
equations.  This  method  has  been  modified  to  include  the  effects  of  an  operat- 
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ing  propeller.  The  details  of  the  basic  method  for  bodies  without  a  propeller 
have  been  provided  by  Chen  and  Patel  (1985)  and  will  not  be  reproduced  here. 
However,  for  completeness,  and  to  aid  in  understanding  the  necessary  modifi¬ 
cations  and  procedures  for  the  incorporation  of  an  operating  propeller,  a 
brief  review  is  given  below. 

The  method  of  Chen  and  Patel  solves  the  transport  equations  for  momentum, 
turbulent  kinetic-energy  k,  and  its  dissipation-rate  £,  in  conjunction  with 
the  continuity  equation,  for  the  mean  velocity  components  (U,V,W),  pressure  p, 
and  turbulence  parameters  k  and  e.  The  continuity  equation  and  the  five  trans¬ 
port  equations  for  (U,V,W,k,£)  are  written  in  cylindrical  polar  coordinates 
(x,r,9)  in  the  physical  domain  as  follows: 
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(x,r,6)  are  the  dimensionless  coordinates  with  (x,r)  normalized  by  a  char¬ 
acteristic  length  L,  and  t  is  the  time  normalized  by  L/Uo.  U,V,W  are,  respec¬ 
tively,  the  longitudinal,  radial,  and  circumferential  components  of  mean  velo¬ 
city  normalized  by  the  characteristic  velocity  UQ;  p  is  the  pressure  norma- 
2 

lized  by  PUq;  Re  =  U^L/v  is  the  Reynolds  number  defined  in  terms  of  UQ,  L  and 

molecular  kinematic  viscosity  v;  the  barred  quantities  uu,  uv,  etc.  are  the 

2  2 

Reynolds  stresses  normalized  by  Uq;  =  c^k  /e  is  the  eddy  viscosity;  Reff 

is  the  effective  Reynolds  number, 
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in  which  <J>  “  k  for  the  k-equation  (II-5)  and  $  =  £  for  the  e  -  equation  (II- 
6);  and  G  is  the  turbulence  generation  terra, 
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The  terras  jib  -  (fb  ,  0,  fb^)  in  the  x-  and  z-momentura  equations,  respectively, 
are  the  components  of  the  body  force  normalized  by  PU2/L  which  represent  the 
influence  of  the  propeller.  These  will  be  discussed  subsequently  in  Section 
II.  B. 

In  the  two-equation  k-e  turbulence  model  used  in  the  present  study,  each 
Reynolds  stress  is  related  to  the  corresponding  mean  rate  of  strain  by  the 
isotropic  eddy  viscosity  v  as  follows: 
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The  model  constants  are: 


—  9V  2 

-  W  -  vt  (2  -  J  k 


-«5-v  C2  |W  +  2  V 

t  r  30  r ' 


=  0.09,  Cgl  =  1.44,  Cg2  =  1.92,  =  1.0,  and  ag  =  1.3 

In  the  present  approach,  the  unsteady  turbulent  flow  is  treated  by  the 
so-called  quasi-steady  approximation  insofar  as  the  basic  k-e  turbulence  model 
derived  from  steady  flows  is  retained.  In  other  words,  any  direct  interaction 
between  the  turbulence  structure  and  the  unsteadiness  induced  by  the  propel¬ 
ler,  or  its  wake,  is  ignored.  Little  experimental  or  theoretical  information 
is  available  at  the  present  time  to  quantify  this  interaction  for  such  complex 
flows.  Generally,  it  is  believed  that  the  quasi-steady  approximation  is  valid 
if  the  frequency  parameter  8  =  f6/U  (where  f  is  the  imposed  frequency 
and  6  is  a  characteristic  boundary-layer  thickness)  is  small  compared  to  a 
characteristic  turbulence  frequency,  such  as  the  turbulent-burst 
frequency,  which  is  of  the  order  of  .2.  This  has  been  observed  experimentally 
for  simple  boundary  layers  with  periodic  external  flows  (Telionis,  1981;  Cous- 
teix,  1986).  For  the  results  to  be  presented,  the  value  of  the  frequency 
parameter  is  6  ~  .03,  and  therefore,  a  quasi-steady  turbulence  model  appears 
to  be  justified,  at  least  as  a  first  approximation. 

The  governing  equations  ( II— 1 )  through  C II— 9 )  are  transformed  into  numer¬ 
ically-generated,  body-fitted,  nonorthogonal,  curvilinear  coordinates  such 
that  the  computational  domain  forms  a  simple  rectangular  region  with  equal 
grid  spacing.  The  transformation  is  a  partial  one  since  it  involves  the 
coordinates  only  and  not  the  velocity  components  (U,V,W).  Referring  to  figure 
1,  for  the  general  case  of  a  three-dimensional  body  the  specified  boundaries 
of  the  solution  domain  are  the  hull  surface  S^,  the  inlet  plane  S^,  the  exit 
plane  Sg,  the  symmetry  planes  Sw  (wake  centerplane)  and  Sj  (waterplane) ,  and 
the  outer  boundary  SQ.  The  outer  boundary  is  placed  at  a  sufficiently  large 
distance  from  the  body  so  that  uniform-stream  conditions  can  be  assumed  there. 
This  eliminates  the  need  of  both  a  potential-flow  solution  and  a  viscous- 
inviscid  matching  procedure.  In  other  words,  the  solution  domain  encompasses 
regions  I  and  IV  of  figure  1. 
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The  transport  equations  are  discretized  using  a  staggered  grid  and  the 
finite-analytic  method  which  reduces  numerical  diffusion  and  instability  as 
compared  with  more  conventional  methods.  The  velocity-pressure  coupling  is 
accomplished  through  the  use  of  a  two-step  global  pressure-correction  pro¬ 
cedure  derived  from  the  equation  of  continuity.  The  equations  are  written  and 
solved  in  the  unsteady  form.  For  steady  flow  applications  time  simply  serves 
as  an  underrelaxation  iteration  parameter. 

The  partially-parabolic  assumption  is  made  in  the  transformed  coordi¬ 
nates;  that  is,  second-derivatives  of  the  transport  quantities  along  the 
longitudinal  coordinate  are  neglected.  As  a  result  of  this  assumption,  only 
the  pressure  field  p  is  fully  elliptic  and  the  velocity  field  is  elliptic  only 
in  the  transverse  plane.  A  solution  may  be  obtained  by  marching  in  the  down¬ 
stream  direction  if  p  is  known  or  assumed;  however,  multiple  global  iterations 
are  necessary  for  the  influence  of  the  downstream  pressure  to  propagate  up¬ 
stream.  The  boundary  conditions  for  a  body  without  a  propeller  are  as  fol¬ 
lows  : 

(1)  on  the  inlet  plane  (S^),  (U,V,W,k,e)  are  given  from  the  thin  boundary- 
layer  (region  III)  and  inviscid-f low  (region  I)  solutions. 

(2)  on  the  outer  boundary  (SQ),  the  uniform-flow  condition  is  applied, 
i.e.,  U  =  U0,  W  =  k  =  e  =  0. 

(3)  on  the  exit  plane  (Se),  a  zero-pressure-gradient  condition  is  applied, 

3p/9x=0. 

(4)  on  the  wake  centerplane  and  waterplane  symmetry  planes  (Sw  and  S^),  W 

=  o,  (U,v,k,e)  =  0. 

(5)  on  the  body  surface  (S^),  a  two-point  wall-function  approach  is  used 
to  specify  (U,V,W,k,e). 

For  a  body  with  a  propeller  the  only  other  additional  boundary  condition 
required  is  that  on  the  propeller  hub  Sp^:  W  =  coR^  ,  where  is  the  hub 
radius.  The  hub  boundary  condition  is  applied  by  using  the  relative  velocity 
in  the  "wall-function"  approach  used  in  the  basic  method.  Thus,  it  is  assumed 
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that  the  log- law  remains  valid  for  a  flow  with  a  spinning  boundary  provided  it 
is  referenced  to  the  boundary.  This  assumption  is  justified  on  the  basis  of 
the  experiments  of  Lohmann  (1976)  even  for  large  values  of 

B.  Representation  of  the  Propeller  and  Interaction  Procedure. 

As  discussed  in  Stern  et  al.  (1985),  there  are  a  number  of  possible 
approaches  for  representing  the  propeller  of  region  II  in  the  partially 
parabolic  formulation  of  regions  I  and  IV.  In  the  present  approach,  the 
propeller  is  represented  not  by  its  actual  physical  presence  but  rather  by  a 
body-force  distribution.  In  this  approach,  following  Schetz  and  Favin  (1977, 
1979),  and  others,  the  influence  of  the  propeller  is  accounted  for  by  the 
addition  of  body  force  terms  fb_  =  _fb_  (x,r ,9 ,t)  in  the  source  functions  of  the 
momentum  equations  (see  equations  (II-2)  and  (II-4)).  It  should  be  recognized 
that  the  body-force  distribution  is  not  a  boundary  condition  but  rather  a 
discrete  force  field  embedded  in  the  flow  field  which  represents  the  effects 
of  the  propeller.  Modifications  to  the  turbulence-model  equations  for  the 
effects  of  the  body  force  (or  rotation)  have  not  been  made  although,  in  prin¬ 
ciple,  they  can  be  readily  introduced.  The  body-force  distribution  is  ob¬ 
tained  from  the  region  II  calculation  for  specified  propeller  geometry,  oper¬ 
ating  conditions,  and  effective  inflow  velocities.  The  vortex— lattice  lift 
ing-surface  method  developed  by  Kerwin  (Kerwin  and  Lee,  1978)  for  predicting 
propeller  performance  or  more  simple  propeller  models  are  used  for  this  pur¬ 
pose.  Since,  in  general,  the  body-force  distribution  depends  on  the  effective 
inflow  velocities  (i.e.,  the  viscous-flow  solution)  which,  in  turn,  depends  on 
the  body-force  distribution,  the  complete  solution  must  be  obtained  itera¬ 
tively. 

The  axial  (fbx)  and  the  tangential  (fb0)  components  of  the  body  force 
are  related  to  the  propeller  thrust  T  and  torque  Q  as  follows: 


R  0  (r)  x  (r,6) 

pie  s 

T(t)  =  /  /  /  fb  dxrdQdr  (11-10) 

Rh  6te(r)  Vr'6) 
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R  9-  (r)  x  (r,9) 
pie  s 

Q(t)  =  /  /  /  fb~  rdxrdOdr 

\  0  (r)  x  (r ,0 ) 

h  te  P 


(II-ll) 


where  the  integration  is  over  the  volume  of  the  propeller  at  its  instantaneous 
position  (see  the  List  of  Symbols  for  the  definitions  of  the  limits  of  inte¬ 
gration)*  Equations  (11-10)  and  (II-ll)  are  the  most  general  form  appropriate 
for  unsteady-flow  simulations*  For  steady-flow  simulation,  it  is  more  appro¬ 
priate  to  use  a  distribution  over  a  disc  (R^  <_  r  <_  Rp)  of  finite  thickness 
(Ax),  i.e. 


R 

rP 

T  =  2ttAx  j  fb  rdr 


(H-12) 


(H-13) 


Results  have  been  obtained  using  both  prescribed  and  interactive  body- 
force  distributions  (Stern  et  al*  ,  1985  and  1986).  In  the  former  case,  the 
overall  solution  procedure  is  noniterative,  i*e.,  the  partially-parabolic 
solutions  are  carried  out  only  once.  This  is  because  equations  (11-10)  and 
(II-ll),  or  (11-12)  and  (11-13),  are  independent  of  the  effective  inflow.  In 
the  latter  case,  several  iterations  are  required,  on  both  the  partially- 
parabolic  and  the  propeller  solutions,  since  fb_  depends  on  the  effective 
inflow. 

For  the  present  investigation,  only  prescribed  body-force  distributions 
were  used.  The  distributions  are  derived  for  conditions  in  which  the  propel¬ 
ler  thrust  coefficient  CT  (or  KT),  torque  coefficient  Kq,  speed  of  advance  J, 
and  radial  circulation  distribution  T(r)  are  known.  In  this  case, 
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(11-14) 


fb 

x 


CtL  R  G(r) 

T  P _ 

R 

4Ax  /  P  G(r)  rdr 


4K_L  R^G(r) 

fb  =  _ 2 _ 2 - 

0  ?  Rn 

7T r J  Ax  /  PG( r)  rdr 

where  L  is  the  length  of  the  body,  and  G  =  r(r)/2rUQR^. 
were  derived  assuming  a  lightly-loaded  propeller. 


(II-15) 


These  distributions 


III.  FLOW  GEOMETRY  AND  EXPERIMENTAL  INFORMATION 

The  most  basic  geometry  for  the  study  of  propeller-induced  flow  fields  is 
a  simple  propeller-shaft  configuration.  This  is  because  the  underlying  trail- 
ing-edge  flow  in  the  presence  of  a  thin  boundary  layer  is  well  documented,  and 
consequently,  it  is  possible  to  identify  and  examine  in  detail  the  important 
features  of  the  propeller— induced  flow  field.  Also,  the  most  extensive  exper¬ 
imental  data  that  is  available  for  propeller-induced  flow  fields  is  for  this 
configuration.  It  is  for  these  reasons  that  this  geometry  was  selected  for 
the  present  study. 

Detailed  laser-doppler  velocimeter  measurements  of  the  flow  field  around 
propeller-shaft  configurations  have  been  made  by  Min  (1978),  Kobayashi  (1981), 
Moyer  and  Billet  (1982),  Billet  and  Davis  (1983),  Okamura  (1983),  and  Wang 
(1985).  The  latter  reference  is  of  particular  interest  since  measurements 
were  made  for  two  distinctly  different  configurations:  propeller  P4660,  which 
is  a  controllable  pitch  propeller  from  the  DD963  class  destroyer  and  has  zero 
loading  near  the  root;  and  propeller  P4498,  which  is  a  research  propeller  and 
has  finite  loading  near  the  root.  P4660  has  a  larger  hub  ratio  (.3)  than 
P4498  (.2),  and  a  lower  pitch  and  camber  at  sections  near  the  root.  Steady, 
circumferentially-averaged,  velocity  profiles  were  obtained  both  in  the  imme¬ 
diate  downstream  vicinity  of  the  propeller  and  in  the  near  and  intermediate 
wake  regions.  The  most  extensive  steady-flow  measurements  are  for  design 
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conditions.  Some  limited  data  was  also  obtained  for  off-design  (large- 
loading)  conditions  and  to  investigate  the  influence  of  hub  length.  Unsteady, 
periodic,  velocity  measurements  were  also  made  in  a  plane  immediately  down¬ 
stream  of  the  propeller.  Here  again  only  limited  experimental  data  are  avail¬ 
able.  The  focus  of  the  experiments  was  on  the  formation  of  the  hub  vortex  and 
its  relationship  to  the  loading  characteristics  near  the  root.  Also,  numer¬ 
ical  predictions  of  the  propeller  circulation  distribution,  using  a  modified 
version  of  the  Greeley  and  Kerwin  (1982)  lifting-surface  propeller  program  to 
include  hub  effects,  were  shown  to  be  in  good  agreement  with  circulation 
distributions  obtained  from  the  measured  swirl  velocity  component  and  Kelvin’s 
circulation  theorem. 

Herein,  results  are  presented  and  comparisons  are  made  with  the  complete 
set  of  experiments  by  Wang  (1985).  Previously  (Stern  et  al.,  1986),  only  an 
overview  of  the  results  was  given.  Also,  improvements  have  been  made  in  the 
unsteady-flow  calculations  such  that  the  solutions  are  fully  converged  and 
extensions  have  been  made  for  off-design  (large-loading)  calculations.  The 
body-force  distributions  used  in  performing  the  calculations  are  based  on  the 
propeller  circulation  distributions  presented  by  Wang.  First,  in  Sections  IV. A 
and  IV. B,  results  are  presented  for  design  conditions  for  P4660  and  P4498, 
respectively.  Second,  in  Section  IV. C,  the  influence  of  hub  length  for  P4660 
at  the  design  condition  is  investigated.  Third,  in  Section  IV. D,  results  are 
presented  for  P4660  for  off-design  (large-loading)  conditions.  Lastly,  in 
Section  V,  results  are  presented  for  unsteady  flow  for  P4660  at  the  design 
condition. 

The  calculation  conditions  were  selected  to  simulate  the  experimental 
conditions  as  closely  as  possible.  However,  this  was  made  somewhat  difficult 
due  to  the  lack  of  sufficient  information  for  viscous-flow  calculations, 
particularly  with  regard  to  the  upstream  flow  (e.g.,  boundary-layer  thick¬ 
ness).  It  should  be  noted  that  the  present  conditions  differ  somewhat  from 
those  used  previously  by  Stern  et  al.  (1986).  Upon  performing  the  calculations 
for  the  off-design  condition,  an  error  was  discovered  in  Table  6.1  of  Wang 
(1985):  the  heading  for  column  9  should  read  Uq  and  not  Kq.  Wang  (1986)  has 
provided  the  correct  values  for  Kq  and  these  have  been  used  in  the  present 
calculations. 
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In  the  presentation  of  the  results  and  the  discussion  to  follow,  all 
coordinates  are  nondimensionalized  using  the  shaft  length  L,  with  x  =  0  at  the 
shaft  nose,  and  velocities  and  pressure  are  normalized  using  the  free-streara 
velocity  UQ  and  fluid  density.  In  many  of  the  figures  the  radial  coordinate 
is  normalized  by  the  propeller  radius  Rp,  and  denoted  by  Y  (=  r/Rp),  and  the 
axial  coordinate  is  also  normalized  by  Rp  and  referenced  to  the  propeller 
plane. 

IV.  STEADY  FLOW 

The  results  to  be  presented  in  this  section  were  obtained  using  the  rota- 
tionally  symmetric  version  of  the  more  general  three-dimensional  partially— 
parabolic  method  described  in  Section  II  (i.e.,  9/99  =  0  in  equations  (II— 1) 
through  (II-9) ) .  This  is  appropriate  for  steady  (circumferentially-averaged) 
flow  around  an  axisyrametric  body  fitted  with  a  propeller.  For  the  design- 
condition  calculations  the  following  conditions  were  employed: 

for  P4660:  CT  =  .746,  KT  =  .279,  Kq  =  .065,  J  =  .976,  to  =  300  rpm, 

for  P4498:  CT  =  .780,  =  .242,  Kq  =  .052,  J  =  .889,  to  =  900  rpm. 

Based  on  these,  and  assuming  a  shaft  length  L=2ra  and  a  water  temperature  of 

15°C,  .the  shaft-length  Reynolds  numbers  for  P4660  and  P4498  were  determined  as 

Re  =  U  L/v  =  2.7  x  106  and  8  x  106,  respectively, 
o 

The  calculations  were  performed  on  a  Prime  9950  minicomputer  and  took 
about  30  minutes  of  cpu  time  for  the  without-propeller  condition.  The  calcu¬ 
lations  were  begun  with  a  zero  pressure  initial  condition  for  the  pressure 

field.  A  time  underrelaxation  factor  oc^_  =  .3  was  used  on  the  body  and  .03 

<_at  <_  .3  for  the  wake.  A  total  number  of  80  global  sweeps  were  used  in  order 
to  insure  that  the  solution  was  fully  converged.  For  the  with-propeller 
condition  a  modification  of  the  above  procedure  was  required  because  the 
addition  of  the  body-force  terms  in  the  momentum  equations  results  in  a  major 
modification  of  the  flow  field  and  can  induce  instability  during  the  iterative 
solution  procedure.  Therefore,  in  order  to  obtain  converged  solutions  for  the 
with-propeller  condition,  the  calculations  were  performed  in  stages,  starting 
with  the  no— propeller  condition  and  building  up  to  the  final  values  of  Crj,  and 
Kq.  In  this  procedure,  the  entire  pressure  field,  as  well  as  the  inlet  veloc- 
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ity  and  turbulence  parameter  distributions,  were  specified  from  the  previous 
solution  along  with  the  increased  values  of  Kq  and  For  the  design- 
condition  calculations,  five  stages  were  required  for  P4660  and  nine  stages 
were  required  for  P4498.  For  the  off-design  condition,  some  further  modifica¬ 
tions  were  required  as  will  be  discussed  in  Section  IV. D. 

A,  Propeller  P4660. 

A  partial  view  of  the  numerical  grid  used  in  the  calculations  for  P4660 
is  shown  in  figure  2.  Referring  to  figures  1  and  2  for  notation:  the  inlet 
boundary  is  at  xu  =  .55;  the  exit  boundary  Se  is  at  x^  =  16;  the  outer 
boundary  SQ  is  at  rQ  -  .4;  the  first  grid  point  off  the  body  surface  is  lo¬ 
cated  in  the  range  50  y+  80  where  y+  =  U^y/v,  is  the  wall-shear  veloc¬ 
ity,  and  y  is  the  distance  from  the  body;  the  number  of  axial  grid  points  is 
85;  and  the  number  of  radial  grid  points  is  48.  This  is  considered  a  fine 
grid  insofar  as  the  first  grid  point  off  the  body  surface  is  as  close  as  is 
allowed  in  the  wall-functions  approach  and  a  sufficient  number  of  grid  points 
is  located  in  the  viscous  flow  region  in  order  to  resolve  it  accurately. 
Turbulent  flat-plate  boundary-layer  profiles  were  used  to  specify  the  initial 
conditions  with  the  inlet  boundary-layer  thickness  estimated  to  be  6  =  .006. 

Calculations  were  performed  for  conditions  corresponding  to  the  shaft 
with  a  fixed  hub  and  no  propeller,  the  shaft  with  a  spinning  hub  and  no  pro¬ 
peller,  and  the  shaft  with  a  spinning  hub  and  propeller.  For  the  with- 
propeller  condition,  two  different  body-force  distributions  were  used.  These 
are  shown  in  figure  3.  Both  of  these  were  derived  from  the  circulation  dis¬ 
tributions  presented  by  Wang  (1985)  using  the  procedures  described  in  Section 
II. B.  One  of  these  corresponds  to  results  from  the  lifting-surface  propeller 
program  modified  to  include  hub  effects  (PSF2-HUB)  and  the  other  corresponds 
to  the  measured  swirl  velocity  component  and  Kelvin’s  circulation  theorem. 
Note  that  there  are  some  differences  between  the  two  body-force  distributions, 
particularly  near  the  hub. 

Figures  4  and  7  provide  an  overview  of  the  results  without  and  with  the 
propeller,  respectively,  through  the  variation  of  some  properties  in  the 
longitudinal  direction.  Figure  5  shows  the  swirl  velocity  profile  for  the  no¬ 
propeller  condition  at  the  hub  apex  (which  is  located  at  x/Rp  =  1*053). 
Figure  6  shows  some  characteristics  of  the  solutions  in  the  far  wake.  Figure 
8  shows  radial  distributions  of  the  various  quantities  of  interest  at  some 
representative  sections  both  upstream  and  downstream  of  the  propeller. 
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Before  discussing  the  influence  of  the  propeller  on  the  flow  field  it  is 
of  interest  to  examine  the  characteristics  of  the  basic  flow,  i.e.,  without 
the  propeller.  Referring  to  figure  4,  it  is  seen  that  for  the  shaft  with  a 
fixed  hub,  the  flow  characteristics  are  consistent  with  a  trailing-edge  flow 
in  the  presence  of  a  thin  boundary  layer.  The  solution  indicates  a  weak  vis- 
cous-inviscid  interaction  and  the  other  usual  features  associated  with  such  a 
flow.  Although  the  flow  on  the  shaft  with  a  spinning  hub  shows  similar  char¬ 
acteristics,  there  are  some  important  differences.  We  observe  a  decrease  in 
pressure  in  the  vicinity  of  the  tail  such  that  the  pressure  recovery  in  the 
wake  occurs  with  an  adverse  as  opposed  to  a  favorable  gradient 5  a  large  in¬ 
crease  in  the  magnitude  of  the  shear  stress  on  the  spinning  portion  of  the 
hub;  the  generation  of  swirl  (W+0);  and  an  increase  in  turbulence  levels  in 
the  near  and  intermediate  wake  (not  shown  here) .  The  axial  and  radial  veloc¬ 
ities  showed  only  a  minimal  influence  of  spin.  This  is  indicated  by  the 
insignificant  change  in  the  wake  centerline  velocity.  Consistent  with  expec¬ 
tation,  the  influence  of  increasing  the  hub  angular  velocity  w  is  a  larger 
decrease  in  pressure,  a  larger  increase  in  the  magnitude  of  the  shear  stress 
on  the  spinning  portion  of  the  hub,  an  increase  in  Wmax,  and  only  rather  small 
changes  in  UCL.  Figure  5  shows  a  comparison  of  the  swirl  velocity  profiles  at 
the  hub  apex  for  three  values  of  the  rotation  parameter  R  =  wR^/U  (=  .97, 
1.26,  1.8)  with  the  experimental  data.  Note  that  the  R  is  related  to  the 

advance  ratio  J  by  R  =  ifR  /JR  .  The  calculations  show  an  increase  in  both 

h  p 

swirl  magnitude  and  radial  extent  with  increasing  R.  Also,  the  radial  posi¬ 
tion  of  W  increases.  The  experimental  data  are  somewhat  scattered  and  only 
max 

roughly  follow  the  calculated  trends. 

Swirling  turbulent  shear  flows  are  classified  as  weak  or  strong  depending 
on  the  magnitude  of  the  rotation  parameter  R  defined  above,  or  a  swirl  number 
S,  defined  by  the  ratio  of  the  axial  flux  of  the  tangential  and  axial  momenta, 

S  =  G./G  R  (IV-1) 

0  x  p 

00 

_  2 

G0  =  /  (pUW  -  puw)  r  dr 
o 

OO 

G  =  /  (pU  +  Puu  +  (p-p  ))  rdr 
x  o 


i  •  0  • 


where 


13 


The  rotation  parameter  is  appropriate  for  shear  driven  flows  such  as  the  shaft 
with  a  spinning  hub  (no  propeller)  while  the  swirl  number  is  more  appropriate 
for  jets  (with  Rp  =  nozzle  radius)  and  propeller  flows.  For  weak  swirl,  i.e., 
R  (or  S)  «  1,  the  axial  velocity  is  not  appreciably  affected  by  swirl. 

Reynolds  (1962)  has  shown  that  the  asymptotic  forms  for  weakly  swirling  wakes 
and  jets  are  identical  to  the  nonswirling  cases  due  to  the  fact  that  the  swirl 
decreases  more  rapidly  than  the  axial  velocity  or  velocity  defect  (see  the 
asymptotic  laws  summarized  on  figure  6).  This  has  been  confirmed  experimen¬ 
tally  for  swirling  jets  (Rose,  1962)  and  for  wakes  (Liu,  1971).  The  present 
results  are  plotted  in  the  format  of  the  asymptotic  velocity  decay  laws  in 
figure  6,  in  which:  Dg  =  shaft  diameter,  Dp  =  propeller  diameter, 

AUCL  =  Uo_UCL  is  the  centerline  velocity  defect,  and  Au  =  U  -U  is  the 
velocity  excess.  It  is  seen  that  the  axial  velocity  defect  follows  the  x“2/3 
power  and  the  swirl  velocity  decays  as  x_1 ,  as  indicated  by  asymptotic  analy¬ 
sis.  Experimental  studies  indicate  that  the  axial  velocity  defect  (or  excess) 
is  decreased  by  swirl,  presumably  due  to  an  increase  in  entrainment  and  mix¬ 
ing.  It  should  be  recognized  that  such  behavior  is  related  to  the  structure 
of  turbulence;  since,  for  laminar  flow,  the  axial  and  swirl  velocities  develop 
independently.  The  present  results  with  the  fixed  and  spinning  hub  only 
roughly  follow  the  experimental  trend  with  regard  to  the  dependence  of  the 
axial  velocity  defect  (or  excess)  on  the  rotation  parameter  (or  swirl  number). 
This  is  presumably  because  the  swirl  is  relatively  small  and  the  basic  k-e 
turbulence  model  used  here  does  not  explicitly  contain  any  effects  of  rota¬ 
tion. 


Next,  the  influence  of  the  propeller  is  examined.  Referring  to  figure  7, 
it  is  seen  that,  as  expected  from  physical  considerations,  a  major  influence 
of  the  propeller  is  an  increase  in  pressure  across  the  propeller  plane,  which 
is  associated  with  the  propeller  thrust,  and  a  large  decrease  in  pressure 
along  the  wake  centerline,  which  is  due  to  the  large  propeller-induced  swirl 
velocity.  The  distribution  of  the  wall-shear  velocity  U  on  the  body  indi¬ 
cates  an  increase  over  both  the  fixed-  and  spinning-hub  conditions  due  to  the 
propeller.  Also,  a  sharp  decrease  is  seen  at  the  propeller  plane.  The  wake 
centerline  velocity  also  indicates  an  increase  over  the  no-propeller  condi¬ 
tions  and  an  overshoot,  i.e.,  values  greater  than  the  freestream  velocity,  are 
observed.  The  swirl  velocity,  which  is  of  course  zero  in  the  absence  of  a 


14 


propeller  and  for  the  fixed-hub  condition,  attains  quite  large  values  just 
downstream  of  the  propeller,  but  decays  with  downstream  distance.  Finally,  we 
note  that  the  pressure  variations  are  quite  sensitive  to  the  body-force 
distribution.  The  results  for  the  PSF2-HUB  body-force,  which  has  larger 
values  of  fb a  and  fbv  near  the  hub,  indicate  a  larger  reduction  in  pressure 
and  wall-shear  velocity  and  larger  swirl  velocities  in  the  near  wake,  than  the 
results  for  the  experimental  body  force.  It  is  believed  that  the  latter  is 
more  accurate,  since  for  this  application,  PS F 2 -HUB  appears  to  overpredict  the 
circulation  near  the  hub. 

Figure  8  shows  the  radial  distributions  of  the  axial,  radial,  and  tangen¬ 
tial  velocity  components  (U,V,W),  pressure  (p),  and  turbulence  parameters 
(k,e)  at  various  stations  both  upstream  (x/Rp<0)  and  downstream  (x/Rp>0)  of 
the  propeller.  Shown  on  figure  8  are  the  results  for  the  conditions  corre¬ 
sponding  to  the  shaft  with  a  spinning  hub  and  no  propeller,  and  the  shaft  with 
a  spinning  hub  and  propeller,  and  for  both  the  experimental  and  PS F 2 -HUB  body- 
force  distributions.  Recall  that  the  solution  for  the  former  condition  is 
similar  to  the  solution  for  the  shaft  with  a  fixed  hub  except  for  the  previ¬ 
ously  noted  differences.  Note  that  x/Rp  -  .28  (figure  8d)  is  just  downstream 
of  the  propeller  trailing  edge  and  x/Rp  =  1.053  (figure  8m)  is  at  the  hub 
apex.  Comparisons  are  made  with  the  experimental  data  wherever  possible.  The 
experimental  data  is  for  the  mean  velocity  components  (U,V,W)  only;  the  cor¬ 
responding  turbulence  measurements  were  not  made.  We  see  that  the  calcula¬ 
tions  are  in  excellent  agreement  with  the  experimental  data  for  all  three 
velocity  components  (figures  8d“8r)  although  some  differences  are  seen  in  the 
W  profile. 

In  general,  all  the  results  are  consistent  with  the  previous  discussion 
and  with  expectations  from  physical  considerations.  The  calculations  indicate 
a  negligible  influence  of  the  propeller  for  x/Rp  =  -2  (figure  8a).  However, 
the  influence  of  the  propeller  becomes  significant  as  the  propeller  plane  is 
approached  (figure  8b-8c) .  The  U  velocity  component  shows  a  large  increase 
across  the  propeller  plane  and  decays  slowly.  At  x/Rp  =  overshoot  is 

still  large  (figure  8s).  The  V  velocity  component  decreases  due  to  the  pro¬ 
peller-induced  flow  contraction.  Note  that  the  shape  of  the  V  profile  close 
to  the  propeller  is  dominated  by  the  propeller  loading  characteristics  (fig¬ 
ures  8c-8e)  whereas  further  downstream  its  shape  is  more  affected  by  the  shaft 
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trailing  edge  and  transition  into  the  wake  (figures  8f-8p).  The  V  velocity 
component  decays  more  rapidly  than  the  U  component*  At  x/Rp  =  2.052  (figure 
8r)  the  calculations  indicate  that  V  ~  0;  however,  the  experimental  data 
indicate  a  slower  decay  such  that  V  is  small  but  nonzero  at  x/Rp  -  2.052.  The 
propeller-induced  swirl  W  is  negligible  upstream  of  the  propeller  (figures  8a- 
8c)  and  is  large  in  and  downstream  of  the  propeller  plane  (figures  8d-8s). 
The  shape  of  the  W  profile  indicates  that  the  propeller-generated  swirl  and 
the  hub-generated  swirl  initially  develop  independently  but  subsequently  merge 
through  the  action  of  turbulent  diffusion.  This  effect  is  due  to  the  fact 
that  P4660  has  zero  loading  at  the  hub  (see  figure  3).  The  experimental  data 
indicates  basically  no  interaction  and  merging  of  the  propeller  and  spinning- 
hub  swirl  velocities  even  as  far  downstream  as  x/Rp  =  1.252  (figure  8o).  This 
is  rather  surprising  since,  from  physical  considerations,  a  more  rapid  mixing 
would  have  been  expected.  The  W  velocity  component  also  decays  slowly  but 
more  rapidly  than  the  U  component. 

The  pressure  profiles  indicate  a  jump  in  pressure  across  the  propeller 
plane.  This  is  associated  with  the  propeller  thrust.  Just  upstream  of  the 
propeller  (figure  8c)  the  pressure  is  reduced  and  just  downstream  (figure  8d) 
it  is  increased.  The  downstream  increase  decays  and  diffuses  outward  (figures 
8e-8m) .  Downstream  of  the  propeller  near  the  hub  and  wake  centerline  there  is 
a  large  decrease  in  pressure  due  to  the  large  swirl  velocity  (figures  8g- 
8s).  This  effect  can  be  explained  with  reference  to  the  V-momentum  equation 
(II-3).  Near  the  propeller  axis  (small  r)  and  for  large  W  the  radial  pressure 
gradient  3p/3r  is  dominated  by  the  -W^/r  terra.  Note  that  the  decay  of  pres¬ 
sure  in  the  wake  can  be  correlated  with  the  swirl  velocity  profile. 

In  general,  the  influence  of  turbulent  diffusion  is  such  that  the  veloc¬ 
ity  components  are  smoothed  and  diffuse  outward  as  the  flow  evolves  into  the 
wake  (figures  8n-8t).  The  calculations  indicate  that  the  turbulence  par¬ 
ameters  k  and  e  are  greatly  increased  in  the  wake  over  the  values  observed 
without  a  propeller,  especially  near  the  wake  centerline  in  the  near  wake  and 
all  across  the  wake  in  the  intermediate  and  far  wake.  As  mentioned  previ¬ 
ously,  no  experimental  data  are  available  to  verify  the  turbulence  predic¬ 
tions.  The  limited  experimental  data  that  are  available  for  other  configur¬ 
ations  (Ravindranath  and  Lakshminarayana,  1981;  Schetz  and  Jakubowski  1975; 
Jessup  et  al.,  1984;  Kotb  and  Schetz,  1985)  indicate  an  increase  in  turbulence 
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due  to  the  propeller.  This  increase  is  due  both  to  the  larger  mean-velocity 
gradients  in  the  with-propeller  condition  and  to  the  turbulence  production  on 
the  propeller  blades.  However,  we  note  that  only  the  former  effect  is 
included  in  the  present  calculations. 

The  velocity  components  in  the  far  wake  with  the  propeller  are  also  shown 
on  figure  6.  Here  again,  the  correct  asymptotic  forms  are  recovered.  Note 
that  the  swirl  number  calculated  at  x/Rp  =  .28  is  S  =  .12.  The  fact  that  a 
propeller  wake  decays  at  the  same  rate  as  the  nominal  wake  of  a  bare  body  has 
been  confirmed  experimentally  by  Aron  (1960).  However,  we  note  that  the 
analysis  of  Reynolds  (1962)  indicates  that  a  moraentumless  swirling  wake  (the 
practical  circumstance)  decays  at  different  rates  than  the  condition  con 
sidered  here  (figure  6).  This  has  not  been  confirmed  experimentally. 

B.  Propeller  P4498. 

Calculations  were  performed  for  P4498  using  two  different  grids,  i.e.,  a 
fine  and  a  coarse  grid.  As  will  be  discussed  subsequently,  the  reason  for 
using  a  coarse  grid  is  that  it  was  not  possible  to  make  calculations  for  the 
with-propeller  condition  with  the  finer  grid  due  to  the  large  loading  near  the 
root  of  P4498.  The  two  grids  are  shown  in  figures  9  and  10,  respectively. 
Note  that  x/Rp  =  .208  is  just  downstream  of  the  propeller  trailing  edge  and 
x/Rp  =  .84  is  the  hub  apex.  The  fine  grid  was  very  similar  to  that  described 
earlier  for  P4660.  In  the  coarse  grid  the  first  grid  point  off  the  body 
surface  was  located  in  the  range  170  _<  y+  _<  220;  the  number  of  axial  grid 
points  was  70;  and  the  number  of  transverse  grid  points  was  40. 

The  initial  and  flow  conditions  for  the  P4498  calculations  are  similar  to 
those  described  earlier  for  P4660.  This  is  also  true  for  the  derivation  of  the 
body-force  distributions  which  are  shown  in  figure  11.  The  very  different 
characteristics  of  the  body-force  distributions  for  P4660  and  P4498  are  evi¬ 
dent  from  a  comparison  of  figures  3  and  11.  In  particular,  P4498  has  finite 
and  large  values  of  fbx  and  fbQ  at  the  root.  Note  that  for  this  case  the 
differences  between  the  experimental  and  PSF2-HUB  body-force  distributions  are 
minimal.  The  format  for  the  presentation  of  the  results  for  P4498  is  similar 
to  that  described  earlier  for  P4660.  Also,  many  aspects  of  the  two  solutions 
are  similar  and  will  not  be  discussed  again.  The  focus  of  the  discussion  to 
follow  is  on  the  differences  between  the  two  solutions. 
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A  comparison  of  figures  12  and  7  shows  that  the  drop  in  the  pressure  for 
P4498  is  considerably  larger  than  for  P4660.  This  is  the  result  of  the  in¬ 
creased  swirl  velocities  for  P4498,  the  maximum  swirl  velocity  at  the  hub  apex 
being  almost  twice  as  large  than  for  P4660. 

A  comparison  of  figures  13  and  8  shows  the  differences  in  the  detailed 
solutions  for  the  two  propellers.  In  this  case  the  experimental  data  is 
restricted  to  the  swirl  velocity  component  only.  Of  particular  interest  is 
the  development  of  the  swirl  velocity  component  which  is  quite  different  from 
that  shown  previously  for  P4660.  For  P4498  the  propeller  generated  swirl  and 
the  spinning-hub  generated  swirl  do  not  develop  independently,  but  interact 
such  that  a  large  hub  vortex  is  formed  at  the  hub  apex  (figure  13j).  Also,  the 
axial  velocity  component  for  P4498  is  much  fuller  near  the  root  than  for 
P4460.  These  effects  are  a  direct  result  of  the  differences  in  the  loading 
distributions.  Note  that  the  calculated  swirl  number  at  x/Rp  =  .208  is  S  = 
.12. 

As  mentioned  previously,  the  results  with  the  P4498  propeller  were  ob¬ 
tained  using  the  coarser  grid.  Calculations  were  made  using  both  grids  for 
the  spinning-hub  condition.  Results  from  both  these  calculations  are  shown  on 
figure  14.  It  is  seen  that  the  coarse-grid  solution  leads  to  a  large  under 
prediction  of  the  maximum  swirl  velocity  and  to  an  overprediction  of  the 
radial  extent.  Note  that  both  solutions  predict  a  considerably  slower  decay 
I n  the  near  wake  (x-  1.307)  than  the  experiment.  Thus,  if  it  had  been  pos¬ 
sible  to  obtain  a  fine-grid  solution  with  the  propeller  for  P4498,  it  is 
expected  that  an  even  closer  agreement  with  the  experimental  data  would  have 
resulted.  The  difficulties  associated  with  the  fine-grid  solution  for  P4498 
can  be  attributed  to  the  use  of  "wall-f unctions1'  in  the  turbulence  model.  The 
present  results  underscore  the  need  for  further  work  on  near-wall  turbulence 
modeling  in  order  to  resolve  the  details  of  the  solution  in  this  region. 
Lastly  for  P4498,  figure  15  shows  the  influence  of  rotation  parameter  for  the 
spinning-hub  condition.  The  swirl  velocity  profile  at  the  hub  apex  is  com¬ 
pared  for  three  values  of  R  (=  .52,  .73,  .87).  The  calculations  show  similar 
trends  as  those  discussed  previously  for  the  P4660  shaft-hub  configuration 
(see  figure  5).  Also,  just  as  was  the  case  previously,  the  experimental  data 
is  somewhat  scattered  and  only  roughly  follow  the  calculated  trends.  Note 
that  the  radial  extent  of  swirl  is  greater  for  the  P4660  shaft-hub  configur¬ 
ation  than  P4498  due  the  larger  values  of  rotation  parameter  R. 


18 


C .  Influence  of  Hub  Length. 

Calculations  were  also  performed  for  P4660  at  the  design  condition  and 
with  a  medium  and  long  hub.  The  previous  results,  discussed  in  Section  IV. A, 
are  for  a  short  hub.  Partial  views  of  the  grids  used  in  the  calculations  for 
the  medium  and  long  hubs  are  shown  in  figures  16  and  17,  respectively.  These 
are  very  similar  to  those  described  earlier  for  P4660  and  shown  in  figure  2. 
Also,  the  initial  and  flow  conditions,  as  well  as  the  body-force  distribu¬ 
tions,  were  the  same  as  before. 

First,  the  influence  of  hub  length  on  the  solutions  without  the  propeller 
are  examined.  Referring  to  figures  4,  18,  and  19,  which  are  for  the  short, 
medium,  and  long  hubs,  respectively,  it  is  seen  that  for  a  shaft  with  a  fixed 
hub,  the  influence  of  increasing  hub  length  Is  to  reduce  the  viscous-inviscid 
interaction.  As  the  hub  length  is  increased,  the  body  curvature  near  the  tail 
is  reduced,  resulting  in  smaller  variations  in  pressure,  and  a  more  gradual 
reduction  in  associated  with  the  thickening  of  the  boundary  layer.  The 
wake  centerline  velocity  is  not  influenced  by  changes  in  hub  length.  The 
influence  of  increasing  rotation  parameter  for  the  medium  and  long  hubs  is 
identical  to  that  described  previously  for  the  short  hub.  The  influence  of 
increasing  hub  length  is  only  a  slight  reduction  in  the  minimum  pressure,  the 
maximum  U  ,  and  the  maximum  swirl  velocity.  Figures  20  and  21  show  a  compar¬ 
ison  of  the  swirl  velocity  profiles  at  the  hub  apex  for  two  values  of  the 
rotation  parameter  R  (=  .97,  1.8)  and  for  the  medium  and  long  hubs,  respec¬ 
tively.  The  trends  are  indentical  to  those  described  previously  for  the  short 
hub  with  reference  to  figure  5. 

Second,  the  influence  of  hub  length  on  the  solutions  for  the  with- 
propeller  condition  is  examined.  Referring  to  figures  7,  22,  and  24,  it  is 
seen  that  in  this  case,  the  influence  of  increasing  hub  length  is  a  further 
reduction  in  pressure  in  the  near  wake,  a  reduction  in  in  the  near  wake, 
and  a  slight  increase  in  Wmax  in  the  near  wake.  These  effects  are  probably 
due  to  the  increase  in  the  length  of  the  development  of  the  boundary-layer 
flow  downstream  of  the  propeller  for  the  longer  hubs. 

Lastly,  figures  23  and  25  show  the  solution  profiles  at  two  stations 
downstream  of  the  propeller  for  the  medium  and  long  hubs,  respectively:  x/Rp  = 

.28  which  is  just  downstream  of  the  propeller  trailing  edge,  and  x/Rp  =  1-616 
and  2.577  for  the  medium  and  long  hubs,  respectively,  which  corresponds  to  the 
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hub  apex.  These  figures  should  be  compared  with  the  corresponding  ones  for 
the  short  hub  shown  in  figures  8d  and  8ra,  respectively.  The  trends  shown  on 
figures  23  and  25  are  consistent  with  the  previous  discussions.  Note  that  the 
comparison  with  the  experimental  data  for  the  swirl  velocity  profile  is  not 
quite  as  good  as  that  shown  previously  for  the  short  hub.  This  is  due  to  the 
fact  that  the  body-force  distributions  used  for  the  medium  and  long  hubs  are 
the  same  as  those  used  previously  for  the  short  hub  and  are  probably  not 
accurate. 

D.  Off-Design  Conditions. 

Results  are  now  presented  for  P4660  for  off-design  (large-loading)  condi¬ 
tions.  Although  only  limited  experimental  data  are  available  for  verifying 
the  results  and,  as  will  be  discussed  subsequently,  the  accuracy  is  severely 
limited,  such  calculations  are  useful  for  demonstrating  the  robustness  of  the 
computational  approach. 

The  open-water  characteristic  curves  for  P4660  are  shown  in  figure  26 
(Wang,  1986).  The  previous  results  were  for  the  design  condition,  which  is  J 
=  .976.  For  the  off-design  condition  calculations,  the  following  were 
employed : 

CT  =  1.82,  Kt  =  .402,  Kq  =  .0894,  J  =  .750,  u>  -  820, 

CT  =  4.75,  Kt  -  .522,  Kq  =  .1126,  J  =  .529,  a>  =  790. 

The  J  =  .75  condition  is  moderate  loading  whereas  J  =  .529  represents  a  very 
heavily-loaded  condition.  The  body-force  distributions  for  J  =  .529  were  de¬ 
rived  as  discussed  previously  and  are  shown  on  figure  27.  Note  that  PSF2-HUB 
predicted  substantially  larger  values  of  circulation  near  the  propeller  tip 
than  those  obtained  from  the  measurements  which  accounts  for  the  differences 
shown  on  figure  27.  As  discussed  by  Wang  (1985),  the  measured  circulation  may 
be  inaccurate  due  to  the  low  data  density  at  this  high  w.  Thus,  the  body- 
force  distributions  derived  from  the  experiments  may  not  be  accurate.  Addi¬ 
tional  inaccuracies  may  also  result  from  the  use  of  definitions  (11-14)  and 
(11-15)  which  were  derived  under  the  assumption  of  a  lightly-loaded  condi¬ 
tion.  Wang  did  not  provide  circulation  distributions  for  J  =  .75.  The  body- 
force  distributions  shown  for  J  =  .75  were  derived  using  the  measured  circula¬ 
tion  distribution  for  J  =  .529  in  equations  (11-14)  and  (11-15)  along  with  the 
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correct  values  of  CT,  Kq,  and  J.  A  comparison  of  figures  3  and  27  shows  that 
for  reduced  J  the  propeller  loading  is  substantially  increased,  including 
larger  values  near  the  root  for  both  fbg  and  fb^. 

Some  further  modifications  of  the  basic  partially-parabolic  method  from 
those  described  previously  were  required  for  performing  the  off-design  condi¬ 
tion  calculations.  This  is  because  the  specification  of  large  values  of  fb0 
can  lead  to  numerical  instability,  especially  in  the  near-wake  region,  unless 
corrective  measures  are  taken.  In  order  to  obtain  convergent  solutions  the 
following  procedures  were  investigated:  (1)  different  types  of  source  terra 

linearization;  (2)  underrelaxation  in  the  momentum,  pressure-correction,  and 
pressure  equations;  (3)  changing  the  sequence  of  solving  the  momentum  equa¬ 
tions  from  U  then  V  then  W  to  U  then  W  then  V;  and  (4)  using  improved  initial 
conditions,  i.e.,  specifying  the  entire  pressure  (p),  velocity  (U,V,W),  and 
turbulence  parameter  (k,£,v  )  fields  from  the  previous  solution  along  with 
the  increased  values  for  and  Kq,  while  building  up  the  solution  in  stages 
from  lower  to  the  final  values  of  CT  and  Kq.  The  first  two  procedures  allowed 
for  calculations  to  be  made  for  J  =  .75  whereas  the  basic  program  with  the 
previously  described  modifications  failed  for  J  =  .85.  Both  procedures  three 
and  four  allowed  for  further  reduction  to  J  =  .529.  The  two  solutions  showed 
some  differences.  In  particular,  procedure  three  predicted  somewhat  larger 
swirl  velocities  and  lower  wake  centerline  pressure  in  the  near  wake  than 
procedure  four.  16  stages  were  required  to  obtain  the  J  =  .529  solution  using 
procedure  four,  and  17  with  procedure  three.  Only  the  results  from  procedure 
four  will  be  presented.  Also,  most  of  the  results  shown  are  for  the  exper¬ 
imental  body-force  distribution;  since,  the  results  for  the  PSF2-HUB  body- 
force  distributions  were  quite  similar  except  for  the  swirl  velocity  profile 
at  the  outer  propeller  radii.  Both  results  for  the  swirl  velocity  profile  are 
shown  for  comparison. 

A  comparison  of  figures  28  and  7  shows  the  very  large  influence  of  large 
loading  on  the  solution.  For  J  =  .529,  it  is  seen  that:  the  near-wake  cen¬ 
terline  pressure  is  reduced  by  almost  an  order  of  magnitude;  the  wall-shear 
velocity  increases;  the  wake  centerline  velocity  shows  an  overshoot  of  2.4  in 
the  very  near  wake  followed  by  a  rapid  reduction  to  1.6  at  x  =  1.5,  a  slight 
increase  to  1.75  at  x  =  2.0,  and  finally  a  monotonic  decrease  to  the  free- 
stream  value;  and  the  maximum  swirl  velocity  in  the  near  wake  is  increased  by 
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a  factor  of  three.  The  large  reduction  in  wake  centerline  pressure  is  a 
direct  result  of  the  increased  swirl  velocities,  and  is  responsible  for  the 
previously  described  wake  centerline  velocity  behavior  which  can  be  explained 
as  follows:  for  J  =  .976  (figure  7),  the  initial  wake  centerline  velocity  in¬ 
crease  is  due  to  the  effects  of  the  initial  favorable  pressure  gradient  where¬ 
as  the  subsequent  increase  is  due  to  turbulent  diffusion  (the  3uv/3r  term  in 
equation  (II-2))  which  dominates  over  the  effects  of  the  subsequent  adverse 
pressure  gradient;  for  J  -  .75  and  .529,  the  wake  centerline  velocity  is 
affected  both  by  the  initial  favorable  (large  overshoot)  and  subsequent 
adverse  (large  decrease)  pressure  gradients  with  turbulent  diffusion  not 
dominating  until  x  >  1.5.  Note  that  the  swirl  numbers  calculated  at  x/Rp  = 
.28  for  J  ®  .75  and  .529  are  S  =  .21  and  .28,  respectively.  Swirling  jet 
studies  indicate  reversed  flow  for  S  =  .6  (Gupta  et  al.,  1984). 

Figure  29  shows  the  solution  profiles  at  the  propeller  trailing  edge 
(x/Rp  =  .28)  and  the  hub  apex  (x/Rp  =  1.053).  A  comparison  of  figure  29  with 
the  corresponding  ones  in  figure  8  (figures  8d  and  8m)  again  shows  the  very 
large  influence  of  increased  loading  on  the  solution.  In  particular,  we  note 
the  following  features:  the  overshoot  in  the  U  profile  is  increased  by  a 
factor  of  2  (for  J  =  .529)  and  is  much  fuller  near  the  wall  due  to  the  finite 
root  loading;  the  V  profile  minimum  and  W  profile  maximum  are  increased  by  a 
factor  of  3  (for  J  =  .529);  the  predicted  swirl  velocities  are  larger  than  the 
measured  ones;  consistent  with  figure  28  there  is  a  large  reduction  in  pres¬ 
sure;  and  the  turbulence  parameters  show  a  large  increase. 

A  part  of  the  discrepancy  between  the  predicted  and  measured  swirl  velo¬ 
city  can  be  attributed  to  experimental  error.  As  pointed  out  previously, 
there  were  difficulties  in  the  experiments  for  this  large  w  which  may  have 
affected  the  measurement  of  the  swirl  velocity  profile  and  the  subsequent 
calculation  of  circulation.  Results  are  also  shown  on  figure  29  for  the  swirl 
velocity  profile  obtained  using  the  PSF2-HUB  body-force  distribution.  Note 
the  significant  decrease  in  W  for  the  outer  radii  due  to  the  larger  circula¬ 
tion  predicted  by  PSF2-HUB  near  the  tip  than  that  obtained  from  the  measure¬ 
ments  • 
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V.  UNSTEADY  FLOW 


The  calculations  presented  and  discussed  in  Section  IV  were  for  steady 
flow  and  neglect  the  effects  of  the  individual  blades  and  their  rotation.  In 
other  words,  the  previous  results  correspond  to  an  infinite-blade  actuator 
disk  representation  of  the  propeller.  In  this  section  extension  of  the  pre¬ 
vious  procedures  for  unsteady  flow  are  discussed  and  results  are  presented  for 
the  design  condition  for  propeller  P4660.  There  are  two  principal  effects  due 
to  a  rotating  f ini te-bladed  propeller:  (1)  the  fanning  action  of  the  propeller 
induces  an  oscillating  pressure  field  all  around  the  propeller;  and  (2)  a 
complex  periodic  blade-to-blade  flow  field,  including  the  individual  blade 
boundary  layers  and  wakes,  is  created  in  the  near  wake  of  the  propeller.  The 
former  effect  is  responsible  for  propeller-induced  hull  vibration.  The  latter 
effect  is  important  for  the  prediction  of  propeller  noise  and  the  details  of 
the  propeller  slip  stream.  As  will  now  be  shown,  only  the  former  effect  can 
be  simulated  with  a  body-force  distribution  for  unsteady  flow. 

Recall  that  the  basic  viscous-flow  solution  method  described  in  Section 
II  is  applicable  to  unsteady  three-dimensional  flow.  This  most  general  ver¬ 
sion  which  solves  the  complete  equations  ( II— 1 )  through  (II-9)  is  required  for 
unsteady-flow  simulations.  However,  since  this  was  the  first  application  of 
this  method  for  unsteady  flow  some  modifications  were  required.  First,  in 
order  to  obtain  fully  converged  unsteady-flow  solutions  a  true  iteration  par¬ 
ameter  was  introduced  and  the  use  of  time  for  this  purpose  was  removed. 
Second,  the  specifications  of  the  boundary  conditions  in  the  transverse  plane 
required  modification  from  the  previously  used  symmetry  conditions  on  the 
waterplane  and  keel.  If  the  upstream  flow  is  axisymmetric,  as  in  the  present 
application,  then  periodic  boundary  conditions  are  appropriate  at  blade  spac¬ 
ing  intervals. 

In  general,  for  unsteady  flow,  the  body-force  distributions  are  distri¬ 
buted  over  the  volume  of  the  propeller  at  its  instantaneous  position  and  are 
time  dependent  as  indicated  by  equations  (11-10)  and  C I I— 11).  However,  con¬ 
sistent  with  the  axisymmetric  inflow  and  the  periodic  boundary  conditions  used 
in  the  present  application,  the  appropriate  body-force  distributions  are  inde¬ 
pendent  of  time  but  are  imposed  over  the  instantaneous  position  of  the  propel¬ 
ler  blade.  In  other  words,  a  rotating,  constant  strength,  body-force  distri¬ 
bution  is  used.  Calculations  were  performed  using  both  a  coarse  and  a  fine 
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cross-sectional  grid  and  propeller  representation.  It  should  be  recognized 
that  the  propeller  representation  and  the  number  of  time  steps  per  period  are 
dependent  on  the  cross-sectional  grid  spacing.  In  the  case  of  the  coarse 
cross-sectional  grid,  the  propeller  was  represented  by  a  single  radial  body- 
force  distribution  (see  figure  30),  and  only  three  time  steps  per  period  were 
used.  In  the  case  of  the  fine  cross-sectional  grid,  the  propeller  was  repre¬ 
sented  by  four  radial  body-force  distributions  (see  figure  31),  and  six  time 
steps  per  period  were  used.  In  both  cases,  the  radial  body-force  distribu¬ 
tions  were  based  on  the  experimental  circulation  distribution  as  previously 
described  (see  figure  3).  In  the  latter  case,  the  chordwise  loading  was  dis¬ 
tributed  based  on  the  NACA  .8  mean  line. 

The  convergence  of  the  solution  for  unsteady  flow  is  influenced  not  only 
by  the  number  of  time  steps  per  period  but  also  by  the  number  of  global  iter¬ 
ations  per  time  step  used  in  the  partially-parabolic  method,  and  the  total 
number  of  periods  calculated.  Presently,  15  global  iterations  per  time  step 
and  a  total  number  of  12  periods  (i.e.,  37  time  steps)  have  been  used  for  the 
coarse  cross-sectional  grid  and  only  2  periods  for  the  fine  cross-sectional 
grid.  The  coarse  cross-sectional  grid  solution  is  fully  converged.  However, 
the  fine  cross-sectional  grid  solution  is  not  due  to  the  limited  number  of 
periods  calculated.  Most  of  the  results  to  be  presented  are  for  the  coarse 
cross-sectional  grid  since  they  are  sufficient  to  demonstrate  the  character¬ 
istics  of  an  unsteady-flow  simulation.  Some  limited  results  are  presented  to 
show  the  effects  of  the  fine  cross-sectional  grid. 

The  initial  and  flow  conditions  used  for  the  unsteady  calculations  were 
identical  to  those  described  earlier  for  steady  flow  and  for  the  design  condi¬ 
tion.  The  radial  grid  distribution  was  somewhat  coarser  than  that  used  for 
the  steady-flow  calculations  (see  figure  2)  but  similar  to  that  described 
earlier  for  P4498  (see  figure  10). 

Figure  32  shows  the  body  surface  and  wake  centerline  pressure  distribu¬ 
tion,  and  wall-shear,  wake  centerline,  and  maximum  swirl  velocities  obtained 
at  the  32nd  time  step  and  for  all  three  grid  lines  over  one  blade  spacing 
interval  (i.e.,  one  period).  These  are  to  be  compared  with  the  corresponding 
steady-flow  results  shown  in  figure  7.  The  unsteady  solution  shows  distinct 
characteristics  at  each  of  the  grid  lines  over  the  blade  spacing  interval: 
for  grid  line  2,  the  propeller  is  in  the  "right  on"  position,  and  its  influ- 
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ence  is  maximum  ,  i.e.,  the  flow  field  variations  are  larger  than  the  steady- 
flow  solution;  for  grid  line  3,  the  propeller  has  "not  arrived"  yet  and  the 
results  are  closer  to  the  no-propeller  solution;  and  for  grid  line  1,  the 
propeller  has  "just  passed"  and  the  solution  shows  a  residual  influence  of  the 
propeller. 

Figure  33  shows  the  radial  distributions  of  (U,V,W,p ,k ,e )  at  some  repre¬ 
sentative  stations  both  upstream  and  downstream  of  the  propeller.  The  results 
shown  are  for  grid  line  2  and  include  all  time  steps  except  the  first  3  (i.e., 
first  period).  The  corresponding  results  for  steady  flow  are  shown  in  figure 
8.  The  results  in  figure  33  clearly  show  the  influence  of  blade  position  as 
described  in  reference  to  figure  32.  Note  that  the  solution  is  fully  conver¬ 
ged,  i.e.,  although  all  periods  are  shown  except  the  first,  only  three  curves 
appear  on  each  figure.  Of  particular  interest,  is  the  extent  of  the  unsteady- 
flow  effects.  It  is  seen  that  upstream  of  the  propeller  the  unsteady-flow 
effects  are  small.  For  x/Rp  =  -1.,  the  unsteady-flow  effects  cannot  be  dis¬ 
cerned  (figure  not  shown).  As  the  propeller  plane  is  approached  x/Rp  =  (-.5, 
-.25),  only  a  gradual  increase  in  unsteady-flow  effects  occurs.  Very  near  the 
propeller  plane  (x/Rp  =  .04)  and  in  the  immediate  slip  stream  .0  <  x/Rp  <  1., 
the  unsteady-flow  effects  are  large.  For  x/Rp  M«»  the  effects  again  become 
small.  For  x/Rp  =  2.,  the  effects  cannot  be  discerned  (figure  not  shown). 
Thus,  the  present  simulation  indicates  that  the  unsteady-flow  effects  are 

limited  to  -1  <  x/R  <  2.  Although  Wang's  (1985)  unsteady-flow  data  are  not 
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extensive  enough  to  confirm  this,  the  earlier  MIT  experiments  of  Min  (1978) 
and  Kobayashi  (1981)  seem  to  indicate  a  similar  extent  (excluding  the  tip 
vortex) . 

Figure  34  shows  the  time  history  of  the  velocity  components  (U,V,W)  at 
x/Rp=.28  for  two  radial  positions  Y=r/Rp=(  .417,  .833)  and  for  grid  line  2. 
Results  are  shown  for  both  the  coarse  and  the  fine  cross-sectional  grid  solu¬ 
tions.  The  results  for  the  coarse  cross-sectional  grid  correspond  to  the  last 
complete  propeller  revolution  (i.e.,  last  5  periods).  Also  shown  for  compari¬ 
son  are  the  experimental  results.  A  blade  angle  coordinate  9  =  ut  measured 
from  the  generator  line  of  the  key  blade  has  been  used  instead  of  time  for  the 
abscissa.  Referring  to  figures  30  and  31,  9  =0  corresponds  to  the  key  blade 
at  grid  line  1.  It  is  believed  that  the  experimental  blade  angle  is  approx¬ 
imately  consistent  with  this;  however,  it  was  not  possible  to  confirm  the 


25 


experimental  blade  angular  coordinate.  Consistent  with  the  steady-flow  pre¬ 
dictions,  the  mean  values  are  well  predicted;  however,  there  are  large  dif¬ 
ferences  in  both  the  magnitude  and  phase  of  the  unsteady  component.  This  is 
especially  true  for  the  radial  V  component.  Note  the  differences  between  the 
coarse  and  fine  cross-sectional  grid  solutions.  The  fine  cross-sectional  grid 
solution  is  clearly  a  more  detailed  representation  and  is  somewhat  closer  to 
the  experimental  trends.  The  trends  shown  by  the  calculations  have  already 
been  described  with  reference  to  figure  32.  The  rotating  body-force  repre¬ 
sentation  of  the  propeller  simulates  the  thrust  and  torque  characteristics  of 
the  rotating  blades  but  is  unable  to  capture  all  of  the  details  of  the  complex 
blade-to-blade  flow  seen  in  the  experimental  data  just  downstream  of  the 
propeller  trailing  edge.  The  experimental  data  reveal  many  interesting  fea¬ 
tures  of  the  blade-to-blade  flow,  including  radial  variations  resulting  from 
the  effects  of  the  propeller  loading,  blade  wakes,  and  tip  vortices.  In  order 
to  predict  such  effects  a  more  detailed  representation  of  the  propeller  is 
required. 

These  unsteady-flow  calculations  were  also  performed  on  the  Prime  9950 
minicomputer  and  took  about  2  hours  of  cpu  time  per  period  for  the  coarse 
cross-sectional  grid  solution  and  4.5  hours  per  period  for  the  fine  cross- 
sectional  grid.  Time  relaxation  was  of  course  not  used,  but  a  relaxation 
factor  of  ot  —  *3  was  used  for  the  pressure  in  the  wake. 


VI.  CONCLUDING  REMARKS 

A  recently  developed  method  for  predicting  propeller-hull  interaction  has 
been  validated  by  performing  comparisons  between  computational  results  and 
extensive  available  experimental  data  for  propeller-shaft  configurations.  The 
steady-flow  results  are  in  excellent  agreement  with  the  experimental  data  and 
show  that  the  present  procedures  are  able  to  predict  most  of  the  details  of 
the  mean  flow  field.  The  dependence  of  the  flow  field  on  propeller  loading, 
including  the  formation  of  the  hub  vortex,  and  the  influence  of  hub  length  are 
accurately  simulated.  Also,  the  robustness  of  the  present  approach  has  been 
demonstrated  by  performing  calculations  for  off-design  (large-loading)  condi¬ 
tions,  including  an  extremely  large-loading  case.  The  results  for  the  off- 
design  conditions  are  quite  realistic,  but  show  differences  from  the  exper¬ 
iments.  It  is  believed  that  a  large  part  of  the  descrepancy  can  be  attributed 
to  the  difficulties  in  the  experiments  for  large  propeller  angular  velocities. 
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Based  on  the  results  from  the  present  investigation  and  those  obtained 
previously  for  rotationally-symmetric  flow  for  DTNSRDC  afterbody  1  and  for 
three-dimensional  flow  for  DTNSRDC  3:1  elliptical  body,  it  is  concluded  that 
the  present  procedures  can  accurately  simulate  the  steady  part  of  the  combined 
propeller-hull  flow  field  and  can  be  used  with  some  confidence  for  predictions 
of  thrust  deduction  and  effective  wake,  as  well  as  other  flow  quantities  of 
interest. 

Of  course,  much  more  work  needs  to  be  done,  especially  for  three- 
dimensional  flow,  to  extend  the  method  to  realistic  ship  geometries.  Some  of 
the  issues  that  need  to  be  addressed  are  as  follows:  optimum  grid  generation 
techniques  in  order  to  resolve  the  flow  field  for  complex  geometries  including 
propeller-induced  effects;  improved  turbulence  modeling  procedures,  especially 
in  the  near-wall  region;  and  the  influence  and  interaction  between  the  propel¬ 
ler-hull  flow  field  and  appendages  such  as  the  rudder. 

The  unsteady-flow  results  obtained  with  the  rotating  body-force  represen¬ 
tation  of  the  propeller  simulate  the  fanning  action  of  a  rotating  finite- 
bladed  propeller.  Such  a  simulation  should  be  useful  for  investigations  of 
propeller-induced  hull  vibration.  Note  that  methods  for  predicting  unsteady 
cavitation  (e.g.,  Stern,  1982)  could  be  incorporated  into  the  present  pro¬ 
cedures  without  too  much  difficulty  to  allow  for  the  inclusion  of  this  impor¬ 
tant  influence  on  propeller-induced  hull  vibration.  The  comparison  between 
the  present  results  and  the  experimental  data  just  downstream  of  the  propeller 
trailing  edge  point  out  the  limitations  of  the  present  approach  for  simulating 
the  complex  blade-to-blade  flow.  A  more  detailed  representation  of  the  pro¬ 
peller  than  the  body  force  is  required.  Two  alternative  approaches  need  to  be 
investigated.  The  first,  is  the  velocity-field  interaction  method  as 
described  by  Stern  et  al.  (1985),  which  is  similar  to  the  present  approach,  in 
that  different  methods  are  used  for  the  propeller  and  viscous  flow  and  the 
complete  solution  is  obtained  iteratively  and  interactively.  The  second,  is 
the  development  of  a  method  in  which  the  actual  propeller  is  embedded  in  a 
viscous-flow  method  and  the  no-slip  condition  is  satisfied  on  the  surface  of 
the  rotating  blades.  It  should  be  recognized  that  both  of  these  approaches 
will  involve  substantial  computer  times  and  the  use  of  a  supercomputer.  Also, 
some  critical  issues  that  need  to  be  addressed  in  the  latter  approach  are  grid 
generation  techniques  such  as  the  use  of  rotating  grids  as  well  as  turbulence 
modeling  and  flow  field  resolution  considerations. 
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Validation  of  the  computational  results  is  made  difficult  by  the  lack  of 
detailed  experimental  data  for  combined  propeller-hull  flow  fields.  This  is 
especially  true  for  turbulence  quantities,  the  slipstream  velocity  field,  and 
the  entire  flow  field  for  three-dimensional  bodies.  Of  particular  interest 
are  detailed  measurements  of  both  the  velocity  and  turbulence  profiles 
upstream  and  downstream  of  the  propeller,  including  mean  and  phase-averaged 
values.  Special  attention  should  be  given  to  the  flow  in  the  immediate  vicin¬ 
ity  of  the  propeller  in  order  to  explicate  the  nature  of  the  unsteady  flow  and 
the  evolution  of  the  blade  wakes  in  the  propeller  slipstream.  An  experimental 
project  is  also  underway  at  the  Iowa  Institute  of  Hydraulic  Research  in  coop¬ 
eration  with  the  Universities  of  Osaka  and  Osaka  Prefecture  of  Japan  with  the 
goal  of  obtaining  such  data. 
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Figure  8.  (continued) 
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Figure  8.  (continued) 
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Figure  8.  (continued) 
42 


00  04  0  8  12  \  6  OO  0,4  0  8  12  16 


U  PROFILE  U  PROFILE 


t  e  profile 


(o) 


(p) 


Figure  8.  (continued) 
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Figure  8.  (continued) 
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Figure  11.  Body-force  distributions 
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Figure  13.  (continued) 
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Figure  13.  (continued) 
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Figure  13.  (continued) 


51 


2  Or 

16  - 


PROPELLER  ^91 

X/R  •  o.a 
» 

- HUB  ONLY 

-  £*P  BOOT  FORCE 

—  —  P5P2-MU8  BODY  FORCE 
o  Exp 


Q6  1.2 

W  PROFILE 


Q* 


Q8  1.2 

W  PROFILE 


I  3  Id 

f  HOFltl 


2* 


(i) 


Figure  13.  (continued) 


52 


-L- 1 - 1 _ I _ l_ 


10  20  30  40 


(k) 


PROPELLER  449* 

r  1  1 

X/Rp.2.26 

- HUB  ONLY 

-  EXP  BOOT  FORCE 

. 

0  EXP 

~  J 

J 

- - 1 - _ 

OO 

2  0, — 


04  0  8 

U  PROFILE 

_T - 1 - 1 - 


01  I.Z 

*  PROFILE 


2  0 

i  i — i  i  ~i - 1  i  i — i - r — 

- p 

— r — i - 1 — 1—1 - 1 - 1 - r — i — t — 

— r 

16 

ti 

- 

1.2 

- 

1  z 

Y 

T 

08 

□  8 

- 

0.4 

J 

i 

CL  4 

i 
/ 1 

1 

l  4 

1 

00 

_ _ 

- 1 - 1 - 1 - 1 - Lr-rr  T  T _ i _ i _ t 

l 

_ L 

0.0 

~  [  ,|-|  1 

— i - 1 — 1 - 1  r  i  i _ j.  , 

i 

1  1 

-16  -1.2  -0  8 


08  12  J  6 

E  PROFILE 


(1) 


Figure  13. 
53 


(continued) 


2.0 


- 1 - 1 — - 1 - 1 - 

PROPELLER  4498 

- 1 - 1 - 

5.26 

- HUB  ONLY 

-  E*P  BOOT  FORCE 

- - PSF  2-HUB  BODY  FORCE 

. 

O  EXP 

J 

/ 

1  / 

- 1 — — i  ......  i  "  i  , 

'G*  -OJ  -Oi  ■ Q 1 

v  PROFILE 


(m) 


-2.0  -16  -1.2  -0  0  -0.4  00 


0  !Q  20  30  40  50  60 

K  (X1000)  PROFILE 


Figure  13.  (continued) 


54 


D,Q  CM  as  I?  is  20 

W  PROFILE 


2.0 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r~ — I — T - t 

1.6  - 


1.2 


-2.0  -1.6  -1.2  -08  -0.4  0.0 

P  PROFILE 

20i - 1 - 1 - 1 - 1 - 1 - 1 - < - < - ■ — n - r 


16 


0  10  20  30  40  30  60 


K  (X1000)  PROFILE 


E  PROFILE 


(o) 


W  PROFILE 


00  0  4  0  8  1  2  1.6  2  0  2  4 

E  PROFILE 

(P) 


Figure  13.  (continued) 


55 


u.ts 

Y 

- 1 - 1 - 

_  X/R  -  0.136 

1  1 - 1 - 1 - 

- FINE  GRID 

-  COARSE  GRID 

1 - 

0.4 

- 

0  EXR 

_ 

—  o 

0.0 

-  1  -  -  I 

i - 1 _ i  i 

J _ 

0.0  0.2  0.4  0.6  0.8 

W  PROFILE 


0.8 

Y 

0.4 


0.0 

0.0  0,2  0.4  0.6  0  8 

W  PROFILE 

0.8 
Y 

0.4 

0.0 

0.0  0.2  0.4  0.6  0.8 

W  PROFILE 


Figure  14.  Swirl  velocity  profile  (without  propeller) 
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Figure  16.  P4660  medium  hub  grid  (partial  view) 


0.4  0 


0.35 


0.30 


0.25  - 


r/L 


0.20 


0.15 


0.10 


0.05 


0.00 


- — - 1 — — h 


0.7 


0.0 


0.9 


1.1 
*/L 
- 1 — 


12 


13 


1.4 


1.5 


2  3  4  5  6 

x/R 


Figure  17.  P4660  long  hub  grid  (partial  view) 
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Figure  18.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (without  propeller) 

59 


X/L 


'•0  3.0  5.0  7.0  9.0  11.0  13.0 


X/L 

Figure  19.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (without  propeller) 
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Figure  20.  Swirl  velocity  profile  (without  propeller) 
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Figure  21.  Swirl  velocity  profile  (without  propeller) 
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Figure  22.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (with  propeller) 
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Figure  24.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (with  propeller) 
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Figure  26.  Propeller  P4660  open-water  curves 


Figure  27.  Body-force  distributions 
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Figure  28.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (with  propeller) 
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Figure  29.  Velocity,  pressure, 
parameter  profiles 


and  turbulence 
(with  propeller) 
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Figure  30.  Unsteady-flow  coarse  cross-sectional  grid 

and  propeller  representation  (partial  view) 
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Figure  31.  Unsteady-flow  fine  cross-sectional  grid 

and  propeller  representation  (partial  view) 
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Figure  32.  Body  surface  and  wake  centerline  pressure, 

and  wall-shear,  wake  centerline,  and  maximum 
swirl  velocities  (with  propeller) 
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Figure  33.  (continued) 
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Figure  33.  (continued) 
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